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Cable theory has been developed over the last decades, usually assuming that the extracellular space around 
membranes is a perfect resistor However, extracellular media may display more complex electrical properties 

p^ ■ due to various phenomena, such as polarization, ionic diffusion or capacitive effects, but their impact on cable 

properties is not known. In this paper, we generalize cable theory for membranes embedded in arbitrarily 

^^ ' complex extracellular media. We outline the generalized cable equations, then consider specific cases. The 

^ ^ , simplest case is a resistive medium, in which case the equations recover the traditional cable equations. We show 

i^ ' that for more complex media, for example in the presence of ionic diffusion, the impact on cable properties such 

^«H . as voltage attenuation can be significant. We illustrate this numerically always by comparing the generalized 

cable to the traditional cable. We conclude that the nature of intracellular and extracellular media may have a 

^"— ^ ■ strong influence on cable filtering as well as on the passive integrative properties of neurons. 

PACS numbers: 

d 

•i-H 

^H. Cable theory, initially developed by Rail [25], is one of the most significant contributions of theoretical neuroscience 

and has been extremely useful to explain a large range of phenomena (reviewed in [26]). However, cable theory makes 
a number of assumptions, one of which is that the extracellular space around neurons can be modeled by a resistance, 
K^ ' or in other words, that the medium around neurons is resistive or ohmic. While some measurements seem to confirm 
1^ . this assumption [19], other measurements revealed a marked frequency dependence of the extracellular resistivity [12, 
\^ \ 13], which indicates that the medium is non-resistive. Indirect measurements of the extracellular impedance also show 
"^ ' evidence for deviations from resistivity [5, 9], which could be explained by the influence of ionic diffusion [6]. Despite 
"^ i such evidence for non-resistive media, the possible impact on cable properties has not been evaluated. 

CO ' The eff'ect of non-resistive media can be investigated by integrating this eff'ect in the impedance of the extracellular 

medium, Zg, and in particular, through its frequency dependence. For example, it can be shown that Z^ ~ 1/w for 
capacitive effects or electric polarization [2], Z,, ~ 1/ V^ for ionic diffusion (also called the "Warburg impedance" [4]), 

^^ I while Zj, would be constant for a perfectly resistive medium. To integrate such effects in a given formalism, such as the 

H ' genesis of extracellular potentials, our approach has been to integrate a aeneral frequency-dependent function Zg( w) in the 
C^ ■ 
- - - formalism, and then consider specific cases [4, 6]. 

In the present paper, we follow this approach and generalize cable equations for media with arbitrarily complex 
frequency-dependent impedance. With numerical simulations, we consider specific cases such as resistive media, ionic 
diffusion, capacitive media, etc. We evaluate a number of possible consequences on the variation of the membrane poten- 
tial along the cable, and how such effects could be measured experimentally. 



II. METHODS 

All simulations were done using MATLAB. To simulate the cable structure of the models, a classic compartmental 
model strategy was used for simulations (see Fig. IIIA2F), but was different from the one used in common simulator 
programs such as NEURON (Hines and Carnevale, [15]). Each cylindric compartment is connected to intracellular and 
extracellular resistances or impedances, and these are normally used to solve the cable equations. In the present paper, we 
used another, equivalent method which consists of defining an auxiliary impedance, given by Za = -^ where V,,, and ;,- 
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are respectively the transmembrane potential and the axial current per unit length at the point where Z„ is connected (see 
Fig. II). This auxiliary impedance allows to take into account the influence of other compartments, including the soma, 
over the axial current and transmembrane potential. It is mathematically equivalent to consider the continuity conditions 
on axial current and transmembrane potential. 

The electric and geometric parameters are considered constant in each compartment, but are allowed to vary between 
compartments. In these conditions, y„, and /, are solution of partial differential equations (cable equations) and thus 
depend on spatial coordinates. 
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FIG. 1 : Convention used to calculate the input impedance and transfer function. A cable segment of length / is represented, with an 
impedance Z^ in series, at the end of the cable. This "auxiliary impedance" Z^ takes into account the influence of the other compartments 
on the axial current ;, and transmembrane potential V,„ in a compact form. Za = ^j4^ where !,(/) is the current per unit length and V,„{1) 
is the transmembrane voltage at coordinate x = I. 

The cable equations simulated in this article are generalized to allow one to include media with complex electrical 
properties. We have designed a MATLAB code that simulates such generalized cable structures, using different types of 
linear density of complex impedances ([Q/m]) and specific impedances ([Q.m]) in each compartment. See Results for 
details of this method. 

All computations were made in Fourier space. We have applied the theory to four different types of media to evidence 
their effect on the spatial and frequency profile of the membrane potential. These models are called SC, FC, FO and NIC, 
respectively (see Table I). The SC model is the "standard model" as defined by Tuckwell [28]; the FC model corresponds 
to a model similar to the standard model (based on a closed circuit), but the cytoplasm and extracellular media impedances 
can be frequency-dependent. The FO type model is the same, with an open circuit (no return current). The NIC model 
includes a non-ideal capacitance similar to a previous study [3]. See Results for details of these models. 

All numerical simulations were made using a "continuous ball-and-stick" model, consisting of a single cylindric com- 
partment, described as a continuum (see Results), and a spherical soma. The dendritic compartment has a radius of 2 /xm 
and a membrane time constant of 5 ms, which corresponds to typical values of in vivo conditions. The has a radius of 
7.5 nm and the specific capacitance was of 0.01 F/rn^. These parameters represent typical values used in a number of 
previous studies [16, 17, 26, 28]. 



III. RESULTS 

We start by generaUzing the cable equations for membranes embedded within extracellular media of arbitrarily complex 
electrical properties. Next, we consider a few specific cases and numerical simulations. 



A. Generalized cable equations 

In this section, we redefine the cable equations taking into account the presence of complex and/or heterogeneous prop- 
erties of extracellular and intracellular media. Because electrically complex or heterogeneous media can display charge 
accumulation, one cannot apply the usual (free-charge) current conservation law. One needs to use a more general con- 
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Types 


Model 


„('") 


A^ 


.,2 _ l + 'tiJT,,, 


SC 


Standard cable 
(closed-circuit) 


Jtfl) Ijnl'e 




(r,+/v)(l+/(LJr,„) 


^^ (7-,+re){l+f(^r,„) 


FC 


Frequency-dependent cable 
(closed-circuit) 


Jm) r„z. 


Zi+Ze 


I'm 


^^ (s,+Ze)(l+f(^r,„) 


FO 


Frequency-dependent cable 
(open-circuit) 


J"i) 


{III) 

i2|L[l + k_(l + /^T„,)] 


Zi{l+i<^T„i) 


Jill) 
rii,[\+~-^([+i(uTn,)\ 


NIC 


Non-ideal cable 
(closed-circuit) 


„(/") ^ niiriiiT^ 


r,„ r{\+ili)T,„){\+iLdTM)-i 

Zj '- l+/(L)(r,i,+r^) J 




*■<' [l+/w(r,„+Tj,)][I+;(ur,„] 



TABLE I: Summary of dendritic cable types and parameters. The table gives the parameters z^" , A' and /c,i for different model types. 
The standard model (SC) is the cable model as given by Rail, Koch and Tuckwell [17, 25, 28]. The "frequency-dependent model" 
(FC) correspond to a standard cable (closed circuit), but where the parameters Zi and Ze are allowed to be frequency dependent. In 
the "frequency-dependent open-circuit model" (FO), the current in the extracellular medium is "perpendicular" to the membrane (see 
Fig. 6). The "non-ideal cable" model (NIC) is similar to the standard model, but the capacitance of the membrane is non-ideal, as 
developed previously [3]. z, (see Eq. 14) and Ze are respectively the impedance per unit length of the cytoplasm and of the extracellular 
medium, respectively, for FC type models. We write r, and r^ when the parameters Zs do not depend on frequency (SC type model). 
The parameter z^'"' (see Eq. 22) is used in FO type models. 

servation law based on the generalized current. In Section III A 1 below, we derive this generalized current conservation 
law, while in Section III A 2, we use this generalized conservation law to derive the generalized cable equations. 



1. Generalized current conservation law in heterogeneous media 

In this section, central to our theory, we show that the free-charge current conservation law {j f) does not apply to 
systems with complex electrical properties. Another, more general, conservation law must be used, the generalized 
current conservation law. We derive here the conservation law for the membrane current in arbitrarily complex media, 
starting from first principles. 

Maxwell theory of electromagnetism postulates that the following relation is always valid for any medium: 



VxH 



;^+ 



dt ' 



(1) 



where H is the magnetic field, and j ^ is the current density of free charges, and ^ is the displacement current density. 

Now consider a domain D delimited by a closed surface dD. If we assume that the medium and the field are sufficiently 
regular, then the divergence theorem applies in D, and we have: 



(rh'^xHhdS = fffw(VxH)dvsO 



(2) 



because the following equality always applies: V ■ (V x /?) s [30] 
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We now define the generalized current density j ^ as: 



^=;^+^-;^+;^ o) 



td 



where y is the displacement current density. 

It is important to note that the term ^ = e,,^ is different from zero, even in the vacuum (assuming that the electric 
field varies in time). 

From Eqs. 1 and 2, we have the following identity : 

(ffij^ hdS = fff'!J j^ dv^O, (4) 

which is valid for an arbitrary domain D. 

In the following, we will use the notation / ^ for the generalized current, ; f for the current due to free charges, and / '' 
for the displacement current, across an arbitrary surface S, according to the following definitions: 

' i9 'g ffJ'^-ndS 
s 

if drf jfjf.hdS (5) 

s 

i^'i' ff§-ndS 



Within these definitions, we can write that the generalized current / ^ is conserved at every time and independently of the 
nature of the medium. At every time, the inward current entering a given domain D is always equal to the outward current 
exiting that domain, independently of the homogeneous or heterogeneous nature of the medium. It is also independent of 
the fact that there may be charge accumulation in some elements of volume, because Eq. 4 always applies. 

Note that this generalized current conservation law does not express anything new on a physical point of view, but is the 
charge conservation law expressed as a function of currents . Indeed, taking into account Maxwell-Gauss law ( V ■ Z3 - p^), 
the definition of j ^ (Eq. 3) and the identity given by Eq. 4, we obtain the differential charge conservation law: 

v.^.v.;/.v.^ = v.;/.^.v.;/.|^o (6) 

ot ot at 

In the case of an electric field in a homogeneous and locally neutral medium, we have ^ ' j ^ - --^ =0 because 
there cannot be charge accumulation anywhere. Because the relation V ■ y ^ = applies to any type of medium, we also 
have V • (^) = 0. Thus, in a homogeneous locally-neutral medium, we have two independent current conservation laws: 
one law applies to the free-charge current / ^ and another one applies to the displacement current / '^ . Note that in a 
homogeneous medium / '' is not necessarily negligible, but the application of the current conservation law on / ^ can be 
done independently of the existence of / '' because the two laws are independent. 

This is the framework assumed in the standard cable theory, in which the extracellular medium is resistive and homo- 
geneous, the displacement current i '' is negligible, and there cannot be charge accumulation inside the dendrites nor in 
the extracellular medium. We will see below that these assumptions do not hold for complex extracellular media. If the 
medium is heterogeneous, then charge accumulation will necessarily appear in the presence of an applied electric field. 
Capacitive effects is an example of such charge accumulation. In such a case, the two current conservation laws on / ^ and 
i '' do not apply to every region of space (see Appendix A). However, the generahzed current conservation on / ^ is still 
valid in all cases. 

Thus, to derive cable equations in heterogeneous media, one must use the generalized current conservation law, as done 
in the next section. 
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2. Application of the generalized consen'ation law to cable equations 

To start, we consider a small portion of membrane surface and build a domain in the intracellular side, which is limited 
by the interior surface of the membrane, while the other surfaces of the domain are located inside the cytoplasm (see 
Fig. Ill A 2). 



Exterior surface of the membrane 



Charge + 



Charge ■ 




Interior surface of the membrane 



Domain D 

FIG. 2: Definition of a domain D inside the cytoplasm and adjacent to the membrane. Due to conductance variations in the membrane, 
or due to charged currents, the total charge in domain D varies. One cannot consider that the displacement current across the surface of 
domain D is zero, because this would be in contradiction with Maxwell-Gauss law (see Appendix A). Black circles represent negative 
charges on the interior surface of the membrane, as well as in the cytoplasm, while blue circles indicate positive charges at the exterior 
side of the membrane. 



Using such a definition, in resting conditions, the intracellular side has an excess of negative charges, which are adjacent 
to the membrane. In such a state, we can calculate the free charge density in this domain from Maxwell-Gauss law: 



2(0 



= (U)Dh 

dX) 



dS - est , 



(7) 



where dD is the surface of the considered domain[31] 



Now, suppose that a conductance variation occurs in the domain (for example following the opening of an ion channel). 
This will induce a charged current in domain D and therefore, there will be a variation of the total charge included within 
domain !D, which implies a non-zero displacement current / '' across the surface dD surrounding domain D (without this 
current, the system would be in contradiction with Maxwell-Gauss law; see Appendix A). 



In such conditions, we have: 



-f = #f-"-» 



(8) 



OD 



Can we neglect this cuiTent to study the variations of the membrane potential along the cable? Because it is difficult to 
give a rigorous answer to this question [10, 27], in particular when / '' is non-zero, we consider the generalized current / * 
because this current is conserved independently of / '' (see previous section). This will allow us to treat cable equations 
without making any hypothesis about charge accumulation inside or outside of the cable. 

Moreover, to stay as general as possible, we include a frequency and space dependence of the electric parameters, 
which will allow us to simulate the effect of media of different electric properties, such as capacitive or di ff'usive [2, 4, 6] . 
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In this context, the Hnking equations must be expressed in their most general form [6]: 

D(x, t) — J ^ Siix, t - t)E(x, t) dr 

< 

J fat) = xr[<(f,f-T)£(f,T) 



(9) 



According to this scheme, the generaUzed current density j'-' inside the cytoplasm obeys: 



«y — CX 



[cr'iix, t - t)E{x, t) + Si(x, t - t)—{x, t)] dr 



(10) 



where cr^ix, f) is the intracellular electric conductivity function and e,(x, t) is the intracellular electric permittivity function. 
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FIG. 3: Compartments and equivalent electrical circuits of the membrane and cable segments. A and C depict different configurations in 
a cable of constant diameter, with their respective equivalent electrical circuits shown in B and D. E is the equivalent electrical circuit of 
a membrane compartment of the cable, and F is the equivalent circuit obtained for three compartments. V, is the intracellular potential 
relative to the reference, Ve is the extracellular potential relative to the same reference, z,- is the cytoplasm impedance, r,„/dx and z'e'^/dx 
are respectively the impedances of ion channels and the input impedance of the extracellular medium as seen by the transmembrane 
current, ig is the output current of a cable element in the extracellular medium, and ;, is the axial current. The membrane potential V„, 
equals V, - V,, and may vary according to the position Xj. 



The first term in the integral accounts for energy dissipation phenomena, such as calorific dissipation (Ohm's differential 
law) and diffiision phenomena. The second term represents the effect of charge density variations in the volume elements. 
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In Fourier frequency space, Eq. 8 becomes algebraic. 

J*.*(x, w) = [cr'-{x, ctj) + icoBi{x, cij)]E{x, to) (11) 

Moreover, we have V x £ = 0, which implies E - -W because electromagnetic induction is negligible in biological 
tissue (in the absence of magnetic stimulation[32]). 

If we now consider a one-dimensional cyUndric cable of constant radius a (Fig. Ill A 2A), the generalized current at a 
position X of the cable can be written as: 

/-IT/ 

i'''(x, f) — i^{x, f) ■ {na h) — —na [crUx, cj) + iaJSiix, w)] ■ ——{x, w) (12) 

ox 

where V, is the intracellular voltage difference with respect to a given reference (which can be far away). In the following 
of the text, we will call "compartment" a cylindric cable with constant radius and with uniform electric parameters (see 
Fig. II). It is important to note that this compartment does not need to be isopotential, and the membrane potential will 
depend on the position on the compartment (see scheme in Fig. II). 

If we assume that the impedance (per unit length) of cytoplasm z, can be expressed as: 

7Ta^[cr:(x, o)) + icosiix, w)J 
then the axial current can be written as: 

if(x,co)^--^(x,u) (14) 

Zi ox 

This expression is similar to the traditional cable equation [16, 26, 28], with the exception that the parameter z, is complex 
(with units of [Q/m])[33]. In addition, the transmembrane current /,f, over a cable length dx can be expressed as: 

■xr A w A , T ^ r^ 3V,„(x,f) , crf,,,,^ , dV,„{x,t) {VJx,t)-E,„) 

ijx, t) = i,„ [x, t) dx = Inadx [ C„, + — (V„,(x, t)- E„,)]=dx[ c,„ + ] (15) 

at e at r„ 

where V„, is the transmembrane voltage, E,„ is the resting membrane potential, C„, is the specific membrane capacitance 
(in F/rn^), Cm is the membrane capacitance per unit length (in F/m), cr^, is the electric conductivity (in S /m), l/r,„ is the 
linear density of membrane conductance (in S /m), e is the membrane thickness (in rn) and /,„ is the transmembrane current 
per unit length (in Ajm) (Figs III A 2C-E). Applying the inverse Fourier transform, we obtain: 

i^ix, 0) = /„, (x, 0) dx = j;Vym{y^, 0) - 27r£„5(0)] dx w = 

(16) 
!,-J;(jc, w) = im (x, (jS) dx - [ioJCm + —]V„,{x, ai) dx co ^ 

Note that we assume here that the resting membrane potential E,„ does not depend on time nor on position in the cable. 

Thus, we can see that the Fourier transform of Eq 14 generates a Dirac delta function for null frequency. In the following 
of the text, we consider frequencies diff'erent from zero, because the zero-frequency component of i^ is zero for a signal 
of finite duration, which is always the case in reality. 

In the model above, the expression of the transmembrane current is identical to the generalized membrane current for 
frequencies different from zero. In this case, the generalized current is given by: 

im - A ■ jm - -Ina dx{cr1^ + iajem)'VV - 2na dxicr"^ + iajs,„)— - dx ( h ia>c,„)V,n - dx im - i^ (17) 

t fffi 

where I is the membrane thickness and A is the membrane surface. 

Assuming that the charge variations inside the channels is negligible, then the generalized current conservation law can 
apply to point B in the equivalent scheme (see Fig. Ill A 2 C), and we can write 

i\x + dx, u) — i^{x, CO) — im{x, a>) — i^{x, co) — i^ix, a>) (18) 

It follows that: 

di- ix,co) - —— dx - -i {x,co) - -/„, (x,oj) dx (19) 

ox 
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Using Eqs. 10 and 15, we obtain: 

^ci -—[{cr1(x,a)) + iciJEiix, cij))—-^(x, ai)] = [idJCm H ]V,„(x,cij) 

ox ox r„, 

Applying the partial derivative on the lefthand term, and dividing by jia^ia-'^. + icuiSi), one obtains: 
d^-Vi 1 dio-'. + icoed dVi 1 1 



Note that if the righthand term was zero, then this equation would be identical to the equation describing the electric 



potential outside of the sources [1, 4, 6], because the V operator equals e^-j- in one dimension, in which case the right 



would be equal to V^V,- H — - ■ W,- where y,- = o"^ + iojsi. 

We can simplify Eq. 18 if the cytoplasm is quasi-homogeneous (assuming the scale considered is large compared to 
inhomogeneities due to subcellular organelles), in which case we can consider that the electric parameters of the cytoplasm 
are independent of position x: cr'r(x, u) - cr'^ico) and e,(x, w) = e,((y). This leads to the following expression: 

1 d'^V i 1 d'^V i 1 

-^-l-(-«,w) = J—- : -—-^{x,Cx)) = [lUCm + —]Vm{x,U)) (21) 

Zi OJc na^{cn + iwe,) ox-^ r^ 

in Fourier space. 

If we now assume that the extracellular medium can also be considered as homogeneous (which will be valid at scales 
larger than the typical size of the cellular elements), then we can model the variations of the membrane potential caused 
by the transmembrane current i^^. We can model this effect from the notion of impedance, without making any hypothesis 
on the current field in the extracellular medium. In this case, one can associate to each cable segment dx the specific 
impedance of the extracellular medium, r" , as seen by the transmembrane current, v" has a similar physical meaning 
as r,„, except that it is a complex number in general. In Section IV, we will see that zj," depends on the direction of the 
current field in the extracellular medium. 

Without any loss of generality, we can write in Fourier space: 

y,(x, (jj) = V„,{x, (jj) + zf \w)i„,(^, oj) ■ (22) 



By substituting this last expression in Eq. 19, we obtain 

— [1 + (1 + l(DT„,)] —-r- = [1 + ICOTmWm 

Zi r„ dx^ 

where t,„ = r„,c„. 
Thus, we can write the system in a form similar to the standard cable equation: 

2 i9^y„,(x,a>) . 

A — ^ = K Vi„{x, to) (23) 

where 



(24) 



K^ = 1 + icJT,, 



for a cylindric compartment (see Eq. B8 in Appendix B). It follows that the general solution of this equation in Fourier 
space (w ^ is given by: 

V,„{x, w) = A+(w)e^ + A'(aj)e^^ (25) 

for each cylindric compartment of length / and with constant diameter (see Fig. 11 for a definition of coordinates). For a 
given frequency, we have a second order diff'erential equation with constant coefficients. 



Bedard & Destexhe (April 23,2013) 9 

In general, one can apply Eq. 25 for different cylindric compartments, as in Fig. Ill A2F. In this case, one must adjust 
the different compartments to their specific limit conditions (continuity of V^ and of the current /.* = ~^^ (see Eq. B4 
in Appendix B). 

Note that Eq. 25 is exact for a cylindric compartment of constant diameter. Thus, it is possible to use this property to 
simulate exactly the full cylindric compartment as a continuum with no need of spatial discretization into segments, as 
usually done in numerical simulators. This is only possible if the cylindric compartment has a constant diameter This 
leads to an efficient method to simulate the cable equations. We will refer to this approach as "continuous compartment" 
in the following. 

As mentioned above, the mathematical forms of Eqs. 23 and 25 are identical to that of the standard cable model, but 
with different definitions of A. Thus, we directly see that the nature of the extracellular medium will change the value of 
these parameters, which become frequency dependent. In particular, we see from Eq. 23 that changing these parameters 
will impact on the spatial profile of the variations of V,,,, if the frequency dependence of the ratio ka = jis affected by the 
nature of the medium. Thus, experimental measurement of the spatial variations of V,,, will be able to identify effects of 
the extracellular impedance only if the ratio ka is affected. 

In the next section, we derive expressions to calculate the input impedance Zi„(P) - Zp'^^ and the transfer function of 
the transmembrane voltage Ft(co) - y'^^'"": between two positions in the cable (as a function of the ratio a:^). Later in 
Section IV, we will see that it is necessary to know these quantities to calculate the spatial variation of V,„ and compare 
the standard model with the cable model embedded into complex extracellular media. 



3. Method to solve the generalized cable 

In this section, we present the theoretical expressions which will allow us to calculate the input impedances needed for 
computing the membrane voltage on a cable with varying diameter We consider the input impedance of the membrane, 
as well as the impedance of the extracellular medium, both of which are needed to calculate the spatial profile of the Vm 
in a given cable segment. 

We proceed according to the following steps: 

1. In the previous section, we saw that it is necessary to calculate the ratio Z,„(f ) - 'j"L'^} at the position of the current 
source, to calculate the V,,, produced at that point. 



B 



T 



^out 3 



1 



z , 2 

out 



FIG. 4: Branching cables. The panels A and B respectively represent a branched cable where a dendrite separates into two daughter 
branches, and its equivalent electrical circuit. The equivalent impedance of segment 1 is equal to the input impedances of segments 2 
and 3 (z,,,,, 2 and Zom 3 ) taken in parallel. 

One strategy is, in a first step, to separate the cable into a series of continuous compartments of constant diameter, 
where parameters a (Eq. 12), z, (Eq. 14), r„, (Eq. 15) and zj," (Eq 22) are constant and specific to each compartment. In a 
second step, one calculates the (transmembrane) input impedance Z"^^ - -fk^ at the begin of each compartment by taking 
into account the auxiliary impedance at the end of this compartment, Za - Z"J,' — .'"^ """' = Z"^ (see Fig. II) if there is no 
branching point. At the branching points, the auxiliary impedances are simply equal to the equivalent input impedance 
of « dendritic branches in parallel (where n is the number of "daughter" branches; see Fig. Ill A 3). Thus, because the 
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input impedance at one end is equal to the input impedance of the other compartment connected to this end, one obtains 
a recursive relation (see Eq. B9 in Appendix B): 



yll+l yryn 



Zi„ (^,„Z«+z,„)e ■'"'"+ (^^,Z"-^0 



^'■" ^^'"^ - .,„ (^,, z;;, + 10 e^^..'^ - (^,,z« - 10 ^ ^ 



where 



Thus, we can write 



Zi 



1 + ^(1 + iDT,„) 



Z'l:' ^F[Zl-Zi„,KA„Jn\ 

This leads to the following expression to relate the first to the «'* segment: 



Z«+i = F [...F [F [Zl ■ Zi„K,^,h]; z,„^,„y...; z,„,^i,„/«] (27) 



Note that this algorithm is a generalization of that used to calculate the equivalent resistance for resistances in series. 
Indeed, for resistance in series we have r^,^ - F(...F(ri; ri); r„) where F(ra; ri,) - ra + r\,. The difference between this 
recurrence function and that of Eq. 29 essentially comes from the fact that there is no current leak in a resistance, while 
there is one in a dendritic compartment. 

2. To calculate the profile of V,,, along the cable, one must use the spatial transfer function '^ ,'p^'^} on a continuous 
cylindric compartment of arbitrary length, and calculate the product of the transfer functions between each connected 
compartment. This leads to (see Appendix C and Eq. C3): 

Ft(1, oj; Zl) = -— ,^'^'j:''- . ,, „ (28) 

K^Z^^, cosh{Kxl) + ZiSinhiKj) 



VJP,„0J) _"Y-fV,n(PM,Oj) 



1 1 V,JP:.C0) ^ ' 



VJPuOj) \j V,n(Pi,aj) 



i=\ 

3. To evaluate Zproxtmai we must calculate the first impedance Z]^ which enters the recursive relation 28. This impedance 
corresponds to the impedance of the soma, which is given by: 

Z]„ = Z, + Z,,, (30) 

where Z^ is the soma membrane impedance and Zcs is the cytoplasm impedance inside the soma. This relation is obtained 
under the hypothesis that the soma is isopotential, and the application of the generalized current conservation law implies 
i ^ - 2Tz^ ~ z +'z where V, and Ve are the electric potentials at both sides of the membrane, inside and outside, 
respectively relative to a reference located far-away. 

The impedance of the bilipidic membrane is approximated by a parallel RC circuit where R - R„, is the resistance and 
Tm — RmCm is the membrane time constant. Thus, Z,'„ can be written as: 



tI 



R„ 



Z„ = Z, + Zcs = z 7^— + Zcs (31) 

1 + lUT^ 

Finally, to evaluate Zdistai, we use the "sealed end" boundary condition Z?, = oo. In this condition, we have Z^^ - 
— coth{K,ili) (see Eq. 26). In the case of a single dendritic branch, we can write: 

^distal ^ Zi_ ^^j^^^jy ^ (32) 

ka 
where I is the total length of the cable. 

In the next section, we turn to numerical simulations to investigate passive cable properties in the presence of complex 
media. We consider the most general case, where both the impedance of the extracellular medium and that of cytoplasm 
can be frequency dependent, and determine the respective impact on the spatial profile and frequency content of the 
transmembrane voltage at the level of the proximal and distal ends of the cable. 
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IV. NUMERICAL SIMULATIONS 

The goal of the numerical simulations is here to show how the physical nature of extracellular and intracellular media 
can influence the spatial and frequency profiles of the transmembrane potential. We present simulations of a "continuous 
ball and stick" model, which consists of a continuous cylindric compartment (described by Eq. 25), connected to a 
spherical soma. In this case, the impedance Z,, of the continuous cylindric compartment is the soma impedance (see 
Fig. 11). We do not investigate here the effect of complex dendritic structures, which is left for future studies. Note 
that what we call a "continuous cylindric compartment" actually represents an infinite number of compartments each 
represented by a resistance in series with a parallel RC circ uit (see Fig. Ill A 2F). 

In a first step, we list the different types of models of intracellular and extracellular media that were used. In a second 
step, we present the results of numerical simulations. 



A. Different types of cable models 

We now explain the parameters used for the simulations of the cable presented in Section IV B. 

1. Cable parameters 

Because the cable equation (Eq. 23) is completely determined by the value of k^ for a given frequency, the spatial and 
frequency profiles of the transmembrane voltage are completely determined if the geometry and boundary conditions are 
set. And because k^i is a function of 4 parameters (r,„, t,„, z,, Zg )(Eq. 24) for a given frequency, we have a four-dimensional 
parameter space where the two last parameters (zt, z" ) can be frequency dependent. We will limit our exploration of this 
parameter space by only varying the physical nature of these impedances for realistic values of r^ and t,„, because the 
influence of these parameters has been largely characterized in previous studies [16, 17, 25]. Furthermore, with t,„ and 
o) fixed, the relation ka - -ti^ depends only on A, and thus, hke the classic studies on cable equations, we will use this 
parameter as a main determinant of the cable properties. 

We will explore the generalized cable equations by considering several typical cases: 

a. Standard cable model The first type of model that we will consider is the "standard cable model" (model SC in 
Table I), identical to that considered by Rail, Koch and Tuckwell [17, 25, 28]. In this model, the neuron is a closed system, 
where the inward and outward currents are balanced, forming a closed circuit (see Fig. 6, left). The extracellular current 
flows parallel to the dendrite, as noted previously [28]. This model is equivalent to consider that the field produced by the 
neuron corresponds to an electric dipole configuration. In addition, this model considers that the extracellular medium is 
resistive, or in other words, that the extracellular impedance is a constant. 

In this standard model, the extracellular impedance Zg is either zero (no extracellular resistivity) as in Rail's and 
Koch's formulations [17, 25, 26], or is equal to a constant, which is equivalent to model the extracellular medium by a 
resistance, as in other formulations [28, 29]. Besides its physical non-sense (the extracellular medium considered as a 
supraconductor), using a zero-resistance is usually justified from the fact that the extracellular resistivity is much smaller 
than the membrane impedance. We will see that this justification does not hold if the medium is frequency dependent, 
in which case for some frequency range the extracellular resistivity may be determinant. Thus, to obtain the general 
expression of A and k^ for the standard model, we set Zg - ~ TTTTlffezTn ^'^ '^l- ^^ *-^^^ Table I). 

b. Frequency-dependent cable model The second type of model is an extension of the standard model, where 
the intracellular and extracellular impedances (z,- and te"\ respectively) are allowed to depend on frequency. This 
"frequency-dependent cable model" (model FC in Table I) can account for example for a neuron embedded in capaci- 
tive or diffusive[34] extracellular media, or if the intracellular medium has such properties, or both. In such cases, the 
appropriate frequency-dependent profiles for the impedances must be used. 

In this frequency-dependent model, if t,„ is fixed, the quantity Ze + Zi completely determines the spatial and frequency 
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FIG. 5: Two different cable models for neurons. Left: Closed-circuit model. This is the standard cable model which forms a closed 
system (all inward and outward currents are balanced) and can be described by an equivalent circuit (bottom; shown here for a two- 
compartment model; R^ and Rj are the extracellular and intracellular resistances, respectively). In this model, the current flows parallel 
to the neuron. Right: Open-circuit model. In this more general model, the current is allowed to flow between neighboring neurons, 
or between the neuron and extracellular space, with no necessary condition of local balance (top). In this case, the neuron is modeled 
by an open circuit (bottom), and the current flows "perpendicular" to the membrane. The equivalent circuit is modeled more generally 
with impedances (Z^ extracellular, Z, intracellular) 



profiles of the V,„, and how they deviate from the standard model (see Table I). To explore the effect of the impedances 
Zi + Ze, we consider three typical cases: "resistive", "capacitive" (which is in fact resistive and capacitive in parallel) and 
"diffusive" (which is equivalent to a Warburg type impedance). Such impedances have also been considered in previous 
studies [4, 5, 23]. 

Note that, in order to simulate the standard model, one must necessarily assume that the real part of t" is negative[35], 
which implies that z^" is not a passive impedance per unit length, but is active, and thus requires a source of energy, as 
pointed previously [8, 20]. This point will be further considered in the Discussion. 

c. Open-circuit model In a third type of model, the "Open-circuit" model (FO in Table I), we use a different ap- 
proach. Instead of considering the neuron as a closed system, where all outward currents must return to the neuron, we 
make no hypothesis about the return currents, and allow for example that neighboring neurons exchange currents [36]. In 
this case, one does not need to describe each neuron by a closed circuit, but all neurons are open circuits are are connected 
together (through the extracellular space). Figure 6 shows the current fluxes of the two models, the standard model is a 
closed circuit where the outward currents loop into the inward currents (Fig. 6A), while in the open-circuit model, all cur- 
rents are exchanged with the surrounding medium (Fig. 6B). These two models correspond to different equivalent circuits 



Bedard & Destexhe (April 23, 2013) 



13 



X 10 



CM 

I 




c< 30 



e 



500 

v(Hz) 



1000 




L(m) 



500 

v(Hz) 



1000 




FIG. 6: Spatial and frequency profile of the membrane potential in the cable model with resistive media. A and B respectively show 
the modulus |/c,(| and the phase <I)[/c,i] of /c,i as a function of frequency v for a continuous ball-and-stick model. C. Modulus of the 
transfer function IF7-I as a function of distance L in the dendritic compartment, for frequencies equal to 5, 50, 100 and 150 Hz (see 
corresponding frequencies in A and B). The red curves correspond to a standard cable model (FC, closed-circuit), with r; = 28x 10' Sl/m 
and rg = 18 x 10' Q/m . The black curves correspond to the same model but in an open-circuit configuration (FO model), with 
r, = 28 X 10' £l/m and z" = 0.01 T,„/2naC„, = 0.4 x 10^^ Q..m. The green curves show a non-ideal cable (NIC) model with 
Tm = O.OIt,,,, r, = 28 X 10' fi/m and r, = Q/m. 



(Fig. 6) 

Note that the Open-circuit cable model is practically equivalent to the traditional (closed-circuit) cable model for an 
isolated neuron, if the impedance of the extracellular medium is negligible compared to the membrane impedance. Indeed, 
if z" and Ze tend to 0, then we have (see Table I): 



lim A 



FO 



— = lim A 

Zi z.^Q 



FC 



(33) 



Similar to the frequency-dependent cable model, we will consider the three types of impedances discussed above 
(resistive, capacitive and diffusive) in the simulations of the Open-circuit model. In this case, we separately consider 
the two quantities z, and Ze because these two parameters directly determine the value of A in models of FO type (see 
Table I). Note that in the Open-circuit model, the real part of z~.'" is always positive, so there is no need of any additional 
energy source (see Discussion). 

d. Non-ideal cable model The fourth type of model considered here is the "non-ideal cable model" introduced 
previously [3]. This model postulated that the membrane capacitance is non-ideal, through the use of an additional 
resistance at the arms of the capacitor; this resistance models the fact that there is some inertia time to charge movement 
(or equivalently, a friction). Such a non-ideal capacitance resulted in a shallower frequency scaling, that is a higher 
capacity of the dendritic tree to propagate high-frequency events [3]. Note that in this model, the extracellular medium is 
modeled as a resistance, so in this respect, the non-ideal cable model is equivalent to the standard model. Mathematically, 
the non-ideal cable appears trough the use of zi" (see Table I), which can therefore be viewed as a particular case of 
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FIG. 7: Spatial and frequency profiles of the membrane potential for a model with resistive extracellular medium and diffusive cyto- 
plasm. A and B: modulus \kx\ and phase <I)[ki] of kx as a function of frequency, for a continuous ball-and-stick model. C. Modulus 
of the transfer function IF7-I as a function of distance for different frequencies (same arrangement as Fig. 6). The red curves con^e- 
spond to a model with zero extracellular resistance. The blue curves show models with open-circuit configuration (FO model with 
20 X 10^ fi.m) and diffusive cytoplasm (z, = ,^^^'''_ Cl/m). The black curves show the same model with 

18 X 10' n/m). Note that for the FC model, \Ft\ 



closed-circuit configuration with a resistive extracellular medium (FC model r, 



progressively increases from 5 to 50 Hz, then decreases between 50 and 100 Hz. 

an influence of the extracellular medium on cable properties. Indeed, the non-ideal cable can be shown to be equivalent 
to - or a particular case of - the open-circuit model, where the V,„ corresponds to V, with a far-away reference (see 
Appendix D). We keep this model here for comparison. 

B. Simulation of the different models 

In this section, we present the results of numerical simulations of the models presented in the previous section (see 
Methods). The goal of these simulations is not to be exhaustive in considering all possible combinations of models, but 
present a few typical configurations. The central question is whether the nature of the extracellular medium can have 
determinant impact on cable properties, and for what type of configuration or parameter values does it happen ? 

Analysis of the spatial profiles ofV^ variations 



In this section, we investigate analytically and numerically different particular cases of extracellular and intracellular 
media to determine how the nature of these media affects the spatial and frequency profile of the membrane potential. We 
consider the transfer functions as defined in Table I. The analyses presented here are limited to a ball-and-stick model, 
which allows a better interpretation of the effect of the physical nature of the different media. The effect of complex 
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FIG. 8: Spatial and frequency profiles for a model with resistive cytoplasm and diffusive extracellular medium. Same arrangement 
of panels as for Figs. 6 and 7, but for difl'erent media. The black curves show the behavior of a closed-circuit (FC) type model with 
resistive cytoplasm (r, = 28 x 10' Q/m) and diffusive extracellular space with Warburg impedance {ze = '^^'V Si/m). The red curves 
correspond to a closed-circuit (FC) type model with z, = 28 x 10' il/in and Ze = Cl/m. The blue curves correspond to an open-circuit 



(FO) type model (r, = 28 x 10' Cl/m, zl 



.(in) 



T„, 0.5 _ 20x10-' 

27raCi„ {l+i)^/ill {l+i)y/w 



n.m). 



dendritic tree morphology will be the subject of a future study. To compare the results from the different models, we have 
considered models with identical geometry (see Methods for parameters). 

Resistive models We first considered the "standard model" with resistive intracellular and extracellular media, as well 
as the non-ideal cable model [3] . In Figure 7, we can see that the nature of the cable model (closed-circuit or open-circuit; 
non-ideal) influences the modulus and the phase of k,i, as well as the spatial profile of the transfer function \Ft\- The 
modulus of the transfer function depends more strongly on frequency in the FC model compared to the two other cases 
(Fig.7C), as observed previously [3]. Note that the parameters of the FO and NIC models were chosen such that they are 
equivalent (see Appendix A). 

Capacitive models Next, we considered models where the cytoplasm and extracellular medium are both of capacitive 
(RC-circuit) type. Note that we considered capacitive effects without ionic diffusion, because if both are combined, the 
resulting impedance is of Warburg type. This type of model will be considered next. With purely capacitive media, 
we observed effects very similar to the resistive model shown in Fig. 7, with slight differences only visible for large 
frequencies (greater than about 200 Hz (not shown). The small dimension of organelles (<< 1 //m^) within cells, as well 
as the distance between neighboring cells (~ 30 nm on average) [11, 22] imply that the capacitance values of the media 
are necessary small compared to the membrane capacitance, and thus the purely capacitive effects (without diffusion) are 
likely to be negligible. 

Resistive models with diffusive cytoplasm We next considered models where the extracellular medium was resistive 
as above, but where the intracellular medium (cytoplasm) was diffusive, and described by a Warburg impedance. Figure 7 
shows the spatial and frequency behavior of this model. We can see that the open-circuit (FO) model shows less attenuation 
with distance compared to the closed-circuit (FC) model. Note that these two models give opposite variations when the 
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FIG. 9: Spatial and frequency profiles for fully diifusive cable models. Same arrangement of panels as for Figs. 6-8, but using a 
continuous ball-and-stick model where both cytoplasmic and extracellular impedances are of difi"usive (Warburg) type. The black 
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closed-circuit (FC) type model with z, = "^ ' Q/ra and Ze = £l/m. The blue curves correspond to a closed-circuit (FO) type model 
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then decreases between 50 and 100 Hz. 
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TABLE II: Resonance frequencies of fully diffusive models for different membrane time constants. The resonance frequencies of \k,-i\ 
as a function of the membrane time constant t,„ (see Fig. 10). 



extracellular medium has a zero resistance: in FC type models, \Ft\ attenuates more steeply as a function of distance 
when the extracellular impedance increases, whereas in FO type models, the attenuation becomes less steep. However, 
the spatial profile of \Ft\ also attenuates less with a diffusive cytoplasm compared to a resistive cytoplasm. The latter 
result is expected, because the higher the frequency the more the impedance "short-cuts" the membrane in this case. Note 
that the Warburg impedance used in all diffusive models considered here was applied for frequencies larger than 5 Hz- 

It is interesting to note that in the FC model, a resonance appears around 24 Hz in the modulus of the transfer function 
ka (Fig. 7A). In contrast, the FO model does not display a resonance. 
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FIG. 10: Parameter ka as a function of frequency for fully diffusive models. The red curves correspond to FO and the black curves 
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Resistive cytoplasm with diffusive extracellular medium Next, we considered the opposite configuration as previously, 
namely a resistive model for the cytoplasm, but a dilfusive extracellular medium. Three sets of parameters were c hosen for 
the extracellular space. First, a FO type model with a resistive cytoplasm and a diffusive extracellular medium described 
by a Warburg type impedance (black curve in Fig. 9), and second, a FC type model with similar parameters (blue curve in 
Fig. 8). These two models can be justified if one takes into account the Debye layer at the edge of the membrane [5, 6, 9]). 
The case with a zero extracellular resistance (short-cut) is also shown for comparison (red curve in Fig. 9). The latter 
model represents the same limit case for both FO and FC models, and therefore constitutes the frontier between the two 
families of curves. 

Fully diffusive cable models Next, we have considered the case where both intracellular and extracellular media are 
diffusive. Figure 9 shows the frequency and spatial profiles of the V,,, for such fully diffusive models. Taking the FO 
model with low extracellular impedance {\ze\ = r^ at 1 Hz) leads to large differences with the FC model (Fig. 9, black) 
compared to the FO model (blue) or the FC model with zero extracellular resistance (red). 

We can see that, in FC type models, the larger Ze, the steeper the transfer function attenuates with distance. In contrast, 
in FO type models, larger z'" lead to less attenuation. This paradoxical result can be explained as follows: in FC models, 
Ze plays as similar role as z,, such that for large values of their real part, thermal diffusion attenuates the signal; in FO 
models, z"J plays a similar role as r,„, and large values of |z™| limit the leak membrane current, reducing the attenuation 
with distance. Thus, for large |z"'|, the dendrites become more "democratic" in the sense that the effect of a given input 
will be less dependent on its position on the dendrite. This is only the case for FO models, however. 

As above, the model with zero extracellular resistance represents the same limit case for both FO and FC models, and 
therefore constitutes the frontier between the two models. 

Resonances with diffusive models One interesting finding is that resonances appear in several models using diffusive 
extracellular impedances (Figs. 7 and 9). This type of resonance was studied further in Fig. 10, where one can see that a 
resonance in \ka\ also implies a resonance in \Ft\: The y,„ still attenuates with distance independently of the frequency, so 
that we always have '-jf^ < 0. In addition, Eq. 23 shows that \-g^\ - \ka\\V,„\, so that the quantity \'-g^\ increases when 
\ka\ increases with frequency, which implies that -L^ becomes more negative because this derivative is always negative. 
It follows that \Ft\ attenuates more steeply with distance when \ka\ increases with frequency. Using a similar reasoning. 



Bedard & Destexhe (April 23,2013) 18 

one can show that \Ft\ attenuates less steeply with distance when |a:,i| diminishes with frequency. We conclude that the 
rate of variation of \Fj\ with frequency is always opposed to that of |a:^|. Consequently, the resonance frequency must be 



^ = because we have ^^ < when ^ > and ^ > when ^ 



We also see that the peak frequency of the resonance continuously depends on the membrane time constant (not shown). 
For example, for t„, - 5 ms, the resonance is at about 24 Hz, and for t„, - 20 ms, the resonance is at about 8 Hz (for more 
details see Fig. 10 and Table 11). It is interesting to note that we have observed resonances only in FC type models with 
resistive extracellular media and diffusive cytoplasm (see Fig. 7), but resonances are present in the two types of models 
(FO and FC) when they are fully diffusive. 



V. DISCUSSION 

In this paper, we have introduced a generalization of cable equations to membranes within media with complex or 
heterogeneous electrical properties. We have shown that generalized cable equations can treat a number of problems 
presently not treatable by the traditional cable equations. We have shown that the nature of the extracellular medium has 
a significant influence on fundamental neuronal properties, such as the attenuation with distance and the spectral profile 
of the transmembrane potential. We enumerate below the consequences and predictions of this work, as well as outHne 
directions for future studies. 

A first advance is to properly define cable equations for complex and heterogeneous media. To solve this problem, we 
have introduced the concept of generalized current, and show that the generalized current is conserved in all situations. 
This stands in contrast with the free-charge current, which is conserved only in special cases. For example, if the medium 
is electrically non-homogeneous (with conductive and non-conductive domains), there will be charge accumulation and 
non-conservation of the free-charge current. Thus the traditional cable formalism, which is based on the free-charge 
current, cannot treat this problem. With the generalized current, however, this problem can be treated in a physically 
plausible way, in accordance with Maxwell equations. 

One drawback of generalized cable equations is that they cannot be solved with available neural simulation environ- 
ments, such as NEURON [15], which implements the traditional cable formalism. Consequently, we have developed a 
specific method for the numerical simulation of generalized cables. This method is implementable with traditional simu- 
lation programs, such as MATLAB. Further work would be needed to determine if generalized cable equations could be 
included in neural simulators, as a special case. 

Note that specialized models different from the standard model were introduced relatively recently [21, 23] to include 
aspects which cannot be treated by the standard model. In [23], the cytoplasm was considered as non-resistive but 
capacitive, and was modeled by a RC circuit. It was estimated that this capacitive aspect is important to understand the 
nature of thermal noise in thin dendritic branches. [21] considers the case of the interaction between closely located 
dendritic branches. In this case, the authors study the phenomenon of surface polarization (see also [2]) and evaluate the 
magnitude of the Maxwell- Wagner time of the effective impedance of the extracellular medium, needed to have significant 
influences over the attenuation profile of the V,„. These two studies show that the physical nature of the intracellular or 
extracellular media can have significant influences on cable properties. However, they do represent very particular cases, 
which motivated the present study where we have attempted to consider a broad range of cases, including both intracellul ar 
and extracellular media, as well as ionic diffusion, which was not treated previously. Thus, the present study generalizes 
those prior studies. 

A second main result of this paper was to also generalize the electrical circuit representing neuronal membranes. Inste ad 
of considering the neuron as a closed system, where all outward currents return to the neuron, we have considered the 
more general case which allows current exchange between neighboring neurons, and thus each is represented by an open 
circuit. We have systematically compared open-circuit (FO) models with the traditional closed-circuit (FC) models, and 
found some important differences. FO models have a transfer function that depends much less on frequency and space, 
compared to FC models (see Figs. 7 and 8). 

We also showed that a previously introduced model of non-ideal cable [3] is equivalent to a traditional cable with 
appropriately scaled extracellular resistances (for frequencies smaller than 100 Hz', see Figs. 7 and 8 in [3], as well as the 
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discussion in that paper). 

One of the most important result is the nature of extracellular or intracellular media can have a strong impact on cable 
properties such as the attenuation with distance. We have observed that the nature of the extracellular medium has an 
opposite impact on distance attenuation on FO and FC models. In FO models, larger extracellular impedances lead to less 
attenuation and electrotonically more compact dendrites. The attenuation can be remarkably diminished for fully resistive 
FO models, with only a few percent attenuation (Fig. 9), whereas for FC type models, the opposite was seen, the dendrites 
become more compact for low extracellular impedances. We can say that in these cases, the effect of distal inputs is close 
to that of proximal inputs, and thus the dendrite is more "democratic". It may be that this remarkable property is present in 
some types of neurons to reduce the attenuation of distal inputs, which constitutes another interesting direction to explore 
in future work. 

Another interesting observation is that diffusive extracellular impedances can give rise to resonance frequencies (see 
Figs. 7 and 9), which also appears as a resonance in |a:,)| (Fig. 10). The resonant frequency depends on the membrane 
time constant, and is in the range of 5-40 Hz, which is well within the frequency range of brain oscillations such as theta, 
alpha, beta or gamma rhythms [7]. It is therefore possible that this resonance plays a role in the oscillatory activity of 
populations of neurons. 

Interestingly, we observed that the input impedance of the extracellular medium (zf ) must necessarily be negative 
in the standard model where the medium is resistive. In a closed-circuit configuration, this means that one must neces- 
sarily assume a source of energy, such as an electromotive force. This source of energy can be simulated by a negative 
impedance. This important point was pointed in previous work, where it was called "anomalous impedance" [8, 20]. 
Interestingly, this constraint disappears in the open-circuit configuration. If the current field is open in the extracellular 
medium, then it is not necessary to assume that z" is negative, and there is no need of such a source of energy. 

In conclusion, we think that the generalized cable equations allow one to treat the problem of how neuronal membranes 
behave in complex extracellular and heterogeneous media. Given the possible strong impact of such media as found here, 
future studies should evaluate in more depth whether such media are indeed influential. A possible approach would be to 
find "signatures" of the extracellular medium from the power spectral density of experimentally observable variables, such 
as the membrane potential (for a related approach, see [5]). The direct measurement of the extracellular impedance, at 
present bound to contradictory experimental results [12, 1 3, 19], should give a definite indication whether the generalized 
cable is a necessary approach to accurately model neurons. 



Appendices 

Appendix A: Displacement current, free current and charge accumulation 

In this appendix, we show explicitly that the displacement current / '^ can be used to formally calculate the charge 
variation in a given domain D. Moreover, we show that the displacement current across a closed surface dS which 
surrounds a given domain D is zero when there is no charge variation inside the domain. 

By definition, the density of displacement current (Eq. 5) in frequency space is given by: 

j ''(x, CO) - icije(x, Lij)E(x, w) (Al) 

where a> = 2nf. By applying the divergence on j '' and taking into account Maxwell-Gauss law, we obtain; 

V ■ /"' = ioN ■ (eE) = icopf (A2) 

Thus, we can calculate the amount of free charges in a given domain D from the density of displacement current in 
frequency space. To do this, we have 



Qf(io) = fffp^a co)du^— fffv ■ /''(f, uj)dv= — (fEl'' -hdS = '-^ (A3) 
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where / '' is the displacement current flowing across surface dS. Applying the inverse Fourier transform, we obtain the 
rate of free charge variation in domain D: 

-^{1)^1' it) (A4) 

at 

Therefore, one can say that the charge in the considered volume does not vary of the displacement current across surface 
dD is zero. Finally, because the differential conservation law for free charges implies: 



? ^^^ " UT ^^if^ "^^ " " nr^ ■ ^^^^' '^ ''" =-ffj^-^'^^^ -' ^(^) ' (^5) 



■ df/{x, t) 
d't 

we can then write: 

i f(t) = i \t) + i f{t) = (A6) 

when the surface is closed and when the free charge conservation law applies. 

Thus, the generalized current entering a given closed surface dD is always equal at every time to the generalized current 
exiting dD, even if there is free charge accumulation inside dD. However, this equality does not allow one to deduce if 
there are variations of free charge density inside dD, because the displacement current must necessary be zero across dD 
to have -^ =0 (see Eq. A4). In other words, it is necessary that the displacement current entering dD is equal to the 
displacement current exiting dD to have a constant charge inside dD. Note that in any given circuit, Kirchhoff's current 
law always applies to the generalized current, even if there is charge accumulation inside the circuit, whereas it applies to 
the free charge current only assuming there is no charge accumulation inside the circuit. 

Appendix B: Input impedance of a cable segment in series witli an arbitrarily complex impedance 

In this appendix, we calculate the input impedance of a cable segment of length I when this segment is connected to an 
arbitrary impedance Zg (see Fig. Ill A 3). 

By definition, we have in x = 0: 

, V,n (0, W) 

; /(0, 0)) 

Applying Eq. 25 allows us to directly express V,„ as a function of the cable parameters. We have 

y,„(0, w) = A+(w) e''' + A-{oj) e-'-'', (B2) 

Similarly, applying Eqs. 14, 17 and 22, we obtain: 

1 7^'"' dV 

Zi r,„ dx 

This last expression allows us to express the current at coordinate x = as a function of the cable parameters: 

i'\Q, w) = — [A+(w) e"^'' - A-{oj) e'"''] (B3) 

Zi 



where 



Zi Tzr^ (B4) 



1 + ^(1 + iu>T,„) 
Thus, the expression for the input impedance Z[ is given by: 



= . (^\ . p^kaI + 1 
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We can then evaluate the ratio jr by using the conditions of continuity of the current and of the voltage at point x - I. 
Applying Eqs. 25 and 14 to that point gives: 



V„(l,oj) = A+(w)+A-(w) (a) 
i:\l,aj) = |i[A+(w)-A-(w)] (b) 



Thus, we have 



Z„ = 



V„(l,co) I- + 1 



'" if {I, CO) ^[^-1] 



and we can write 



A+ ^ K^Zg + Zi 

A- KxZa - Zi 



It follows that the input impedance Z\^ is given by: 



'"^ °^ K^ {KAZa+-Zi)e^-''-{KAZa-Zi) 



where 



1 + ^(1 + icjT„i) 



Note that Z'^JZa] — > — when / — » oo, and Z^Za] — » — coth{Kj) when Z,, ^ oo. 



Appendix C: Calculation of the transfer function Ft 



(B6) 



(B7) 



(B8) 



(B9) 



In this appendix, we calculate the transfer function Fjih u; Za) - ^'"1,^^ using the same conditions and conventions as 
for Appendix B. 



Applying Eq. B5a gives: 



y,„(0,w) = A+Me-^'' + A-(w)e-'^-'' 
y„(Z,w) = A+(w) + A-(w) 



Thus, we have 



F7-(^,W;Za) 

Applying Eq. B7 gives the transfer function: 

Fj{l, w; Zn) 
where 



A^(oS)+A-(oS) 
A+(w) e'^j' + A-(w) e 



-,.,,/ 



^iZfl 

A:iZa coshijcxl) + Zisinh(Kxl) 



Zi 



Zi - 



1 + ^(1 + iaiTm) 



(CI) 



(C2) 



(C3) 



Note that Fr(/, w; 0) = and Ft(1, co; oo) 



1 



coshitcj^l) ' 



Bedard & Destexhe (April 23, 2013) 

Appendix D: A new Interpretation of the non-ideal cable 



22 



In this appendix, we show that the non-ideal capacitance model introduced previously [3] is equivalent to an open- 
circuit resistive model if we assume that the circuits A and B in Fig. C are linked by the following transformation: 



^a ~ ^m 



n = 






= fa+rb 

= rb + - 



(Dl) 



(ra+r,,)- ''" 



We show that the V,,, in the non-ideal cable model corresponds to V, in an open-circuit (FO) type resistive model, with a 
reference located far-away. According to circuits A and B in Fig. C, we have: 



j circuit A (r^c ©c,„) || r^ 
1 circuit B (r„ \\ Ca) ffi rb 

It follows that the impedances of circuits A and B are equal if we have: 

Vi2(w) r„ + i(L>r,„r,cC,n r^ + n, + icL)Car„rb 



i{u) 



1 + iojc„(r„, + r,c) 



(D2) 



(D3) 



We see that the ratio -rrrr- is a homographic transform of variable cj. Consequently, Vw we have the following relation 



Vnio') 









i(aj) ~ 1+d oj ~ \+d 01 ^^en the two circuits are equivalent. The only way to guarantee that the equivalence is 
independent of frequency is to assume that the corresponding coefficient of the transformations are equal. We can thus set 
ua - as, bA - bs and cIa - ds- This gives us 3 equations which link the 3 parameters of circuit A to those of circuit B. 
The solution is the transformation law (Eqs. Dl). Thus, on a physical point of view, one cannot distinguish the topology 
of circuits A and B if we would perform external measurements. Moreover, because the functions r„, - fm{ra, rb.rc). 



B 
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> 
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FIG. 11: Equivalence of the electrical circuits of open-circuit and non-ideal cable models. The circuits A and B are equivalent when 



the ratio 






of the voltage difference between points 1 and 2 and the input current between these points is invariant, and when the 



correspondence between the elements of these circuits are independent of frequency. Note that the values of the elements between 
the two circuits are related by a transformation law which is independent of frequency; this equivalence also applies to the temporal 
domain. In other words, according to this equivalence, the two circuits are equivalent when it is impossible to distinguish their topology 
from external measurements. The circuit A conesponds to the non-ideal capacitance model introduced previously [3], while circuit B 
corresponds to a "standard cable model" with a short-cut (zero extracellular resistivity). 
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rsc - fscira, rb-rc), rc,„ - fc„{ra, rb-fc) do not depend on frequency, their equivalence will be also valid for all frequencies. 
We can deduce that the two circuits will behave identically as a function of time . 

It follows that the y„, (between points 1 and 2) in circuit A (non-ideal capacitance) corresponds to the V, relative to a far- 
away reference in circuit B (see Table I). Therefore, a model with non-ideal capacitance and zero extracellular resistance 
should produce a V„, equivalent to the V, of a model with ideal capacitance and resistive extracellular medium. Thus, the 
frequency-scaling behavior of the V„, obtained in a previous non-ideal cable model [3] also applies to the resistive FO 
model, but only if one studies the intracellular potential V,. 
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